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ABSTRACT 



We present images and light curves of the bipolar Planetary Nebula RPJ 053059-683542 that was dis- 
covered in the Reid-Parker AAO/UKST Ha survey of the Large Magellanic Cloud (LMC). The emission 
from this object appears entirely nebular, with the central star apparently obscured by a central band of ab- 
sorption that bisects the nebula. The light curves, which were derived from images from the SuperMACHO 
project at CTIO, showed significant, spatially resolved variability over the period 2002 January through 2005 
December. Remarkably, the emission from the two bright lobes of the nebula vary either independently, or 
similarly but with a phase lag of at least one year. The optical spectra show a low level of nebular excitation, 
and only modest N enrichment. Infrared photometry from the 2MASS and SAGE surveys indicates the 
presence of a significant quantity of dust. The available data imply that the central star has a close binary 
companion, and that the system has undergone some kind of outburst event that caused the nebular emission 
to first brighten and then fade. Further monitoring, high-resolution imaging, and detailed IR polarimetry 
and spectroscopy would uncover the nature of this nebula and the unseen ionizing source. 

Subject headings: Magellanic Clouds — planetary nebulae; general — stars: evolution 



1. Introduction 

Variability in the nuclei of planetary nebulae (PNe) 
has been studied for decades, with renewed inter- 
est in the past few years in order to understand 
the frequency of PNe with binary central stars (e.g. 
Bondll2000l : iMoe fc De Marcoll2006h and the role they 
might play in shaping the host nebula morphology 
and evolution. Variability within the nebulae them- 
selves is less commonly reported, but when present has 
been attributed to a variety of causes. For example 
Feibelman. Aller fc Hvun"3 (1992) attributed changes in 
excitation of the nebulae IC 4997 and NGC 6572 to an 
incr ease in the temperature of th e central stars (CSs); 
and lKato. Nogami fc Babal (|2001l) reported two episodes 
of a deep fading and consequent re-brightening of the 
CS of NGC 2346, which they attributed to dust clouds 
within a PN passing in between the CS and the observer. 



More recently, iDovle et al. 112003) described the remark- 
able morphological evolution and photometric variation 
of M2-9, which they attributed to the interaction of 
ejected matte r between a jet from the C S and the sur- 



round ing PN: lLutz et al.l (|1989 ) and later ICorradi et al 
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(I2OOII) argue from detailed UV and optical spectroscopy 
and IR photometry that IIe2-104 is a PN with a symbi- 
otic (white dwarf plus a Mira variable) central star; and 
Pena et all ( 20041 ) obtained extensive UV and visual- 
band spectra of the complex nebula LMC-SMP 83, 
which they attribute to episodic bipolar ejections from 
a non-degenerate secondary in a close orbit with a white 
dwarf. 

We have initiated a study of photometric variability 
in Large Magellanic Cloud (LMC) PNe in order to make 
use of a very large, complete sample to set limits on the 
fraction of PN central stars that have companions of 
some sort, and to explore the frequency of other pat- 
terns of variability. We have used data from the exten- 
sive SuperMACHO photometric survey at CTIO of the 
LMC which was obtained over a 4. 2 year period from 
2002 through the end of 2005 (see ICarg et"all |2007[) . 
One object, RPJ 053059-683542 (hereafter referred to 
by its common catalog name, RP916), was found to have 
spatially resolved variability, which is a highly unusual 
phenomenon that to date has only one other rough ana- 
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log: the variabilit y in the Galact i c PN M2-9. RP916 
was discovered bv iReid k, ParkeJ ( 20061 ) in the Anglo- 
Australian Observatory (AAO)/UKST deep Ha survey 
of the LMC, who classified it as a "true" planetary neb- 
ula based on its morphology and spectral characteristics. 
It is one of the largest PNe known in the LMC, and has 
a bipolar morphology. Aside from the near uniqueness 
of the phenomenon, this type of variability is highly sug- 
gestive of an interaction of the PN central star with an 
otherwise undetectable binary companion. 

We describe in §2 the extant observations for RP916, 
including the photometric campaign, optical spec- 
troscopy, and broad-band IR photometry from the 
2MASS catalog and from the SAGE survey with the 
Spitzer Space Telescope. We analyze the characteristics 
of the variability in §3, and suggest interpretations that 
are most consistent with the observations in §4. 

2. Observational Data 



Even before iReid &: Parked reported the ini- 

tial discovery, a good deal of ground- and space-based 
data had accumulated on RP 916, including broad-band 
optical and infrared images, and optical spectra. Inter- 
estingly, the bulk of the observations were obtained dur- 
ing the course of surveys of the LMC that were designed 
with purposes other than PN research in mind. 

The primary optical images that were used for this 
study were obtained as a part of the SuperMACHO sur- 
vey of a 23 de g^ regio n centered on th e bar of the LMC 
(seelResTeEal. 2005; Garg et al.ll2007i for details of the 
observing campaign). Eighty- two images in the vicin- 
ity of RP916 were obtained with the Mosaic-II CCD 
camera on the Blanco 4-m telescope at Cerro Tololo 
Inter- American Observatory between 2002 January and 
2005 December. Most images were taken with a broad 
VR filter with a passband that covers 510-740 nm and 
includes the strong nebular emission lines such as Ha 
and [N II] AA6548, 6583, but avoids [O III] A5007. This 
wide bandpass enabled the detection of faint objects 
{mvR ~ 17.5 — 23.5) within the typical 150 s exposures 
of the survey. The CCD plate scale is 0'.'27 pixel"^, 
which samples the PSF very well: the delivered image 
quality varied from roughly O'.'S to 2'.'0 over the course 
of the observing campaigns. An atmospheric dispersion 
corrector was used to minimize the differential atmo- 
spheric dispersion through the broad filter at the typical 
airmass of 1.3 to 1.4. 

The SuperMACH O data were pr o cessed with an 



autom atic pipeline (|Garg et al.l 120071 : iMiknaitis et al 



I2OO7') that corrects for electronic cross-talk between the 
amplifiers, removes the bias level, and applies a fiat-field. 
The images were re-projected to a common geometry, 
with an RMS accuracy of ~ 80 mas. The PSFs are then 
matched, the images are corrected for sky background. 



and placed on a common photometric scale prior to tem- 
plate subtraction. The template image (see §3) is shown 
in the left-hand side of Figure [U and was selected based 
on the good image quality O'.'S) and the excellent 
photometric quality of the night on which this image 
was obtained. 

The discovery images used by iReid k Parked (|2006l ) 
were obtained as part of an AAO/UKST deep photo- 
graphic survey in Ha and in i?— band continuum, and 
were obtained between 1998 and 2000. The 70 A band- 
width of the Ha filter also includes the [N II] emission 
lines at AA6548, 6583. The photometric depth of the 
stacked images is ~ 21.5 in R and Requiv ^ 22 for Ha. 
The digitized images were sampled at 0'.'67 pixel^^, and 
the image resolution is roughly 3'.'5. As described by 



ReidI (|2007^ . the stacked Ha and i?— band images were 
each assigned a specific color in the final image: red for 
continuum, and blue for Ha. The discovery image of 
RP916 is shown in the right-hand side of Figure [H 

Fortuitously, additional images were obtained by one 
of us (A.R.) with Hubble Space Telescope using the 
WFPC2 camera. These images were obtained to study 
a microlensing target, but happened to include RP916 
at the edge of the field of view. In all, four im- 
ages were obtained on 2007 May 17: u9px0601m and 
u9px0602m with the F555W filter, and u9px0603m 
and u9px0604m with the F814W filter. The exposure 
times were all 500 s. We combined the calibrated im- 
ages for each filter to improve the signal-to-noise ratio 
and to remove cosmic-rays. The F814W image is shown 
in Figure m 

Reid fc Parker (I2OO6I) obtained two confirmatory 



spectra of RP916 on 2004 December 14 with the 2dF 
fiber spectrograph on the Anglo-Australian Telescope: 
one with a dispersion of 4.3 A pixel^^ and the other with 
a dispersion of 1.1 A pixel" ^. The higher-dispersion 
spectrum was used to resolve close emission line blends 
and to determine the radial velocity. Since the angular 
size of the nebula exceeds the areal coverage of the fiber, 
the spectra actually correspond to the western lobe of 
RP916. The spectra show moderately low excitation, 
with F([0 III] A5007)/F(H/3) ~ 2.4, relatively strong 
[O I] A6300 and [O II] A3727, weak [O III] A4363, no 
detectable He II emission, and moderate [N II] with 
F(A6548-KA6583)/F(Ha) - 0.5. 

RP916 is a moderately strong infrared source (for a 
PN in t he LMC), based on th e 2MASS J,H,K mag- 
nitudes ( Skrutskie et al.l [2OO6I ) that were obtained on 
2000 February 20 (MJD 51,602). It was also detected 
in the mid-IR SAGE surv ey with Spitzer Space Tele- 
scope ( Meixner et al.|[2006l ) in the 3.6/i, 5.8/i, and 8.0/u 
bands with IRAC, and the 24^ band with MIPS (the 
70/i and 160/i band catalog data have yet to be pub- 
lished). The first epoch SAGE observations were ob- 
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tained during 2005 July 15-26, around MJD 53,570. 
The IR brightnesses are given in Table [T] 

3. Analysis 

All of the optical images show that RP916 has a "but- 
terfly" bipolar morphology ( Balick fc Frank|[2002 ). with 
a dark lane (0'.'46, or 0.11 pc, in width) bisecting the 
pinched main lobes of emission. The extent of the neb- 
ula is 7'.'0 X 3'.'0 (1.72 x 0.74 pc at the distance of the 
LMC) as measured from a contour at 10% of the peak 
nebular brightness in the F814W image; the extent is 
roughly twice as large when measured just above the 
sky level, with a substantial extension of faint nebular 
material to the ENE of the geometric center of the neb- 
ula. The Ha + R color image shows that the emission 
is almost entirely nebular: what little continuum is evi- 
dent in the i?— band image undoubtedly originates from 
a combination of nebular continuum, emission lines in 
the R bandpass, and the few very faint field stars that 
are seen in the HST image. Interestingly, a faint star is 
evident in the F814W image, within the waist of RP916, 
offset by ^ 0'.'23 from the symmetry axis defined by the 
optical emission. The star is barely visible in the F555W 
image however, indicating that it is intrinsically very 
red, or suffers heavy extinction from dust within the 
waist of the nebula. It is not clear whether this is the 
central star of the nebula. The more detailed HST im- 
ages show clumps of emission, with the brightest knots 
or lobes extending almost symmetrically to the east and 
west but do not intersect at the geometric center of 
the nebula. The appearance of RP916 in the 2MASS 
images that of an extended point source, with no sign 
of bi-lobed structure; the object is not resolved in the 
Spitzer IRAC o r MIPS images. Although the confirma- 
tory spectra of Reid fc Parkeri ( 20061 ) are neither deep 
nor well enoug h calibrated to support a detailed abun- 
dance analysis . iReidl (j2007l ) derived an approximate den- 
sity for RP916 of ~ 400 cm-^ from the [S II] A6716, 6731 
line ratio. The amount of extinction, while only roughly 
determined, appears to be higher than average (c ~ 0.9) 
for an LMC PN; given the strong IR emission, some of 
the extinction may be intrinsic to the nebula. Finally, 
the modest ratio of F([N II])/F(IIq;) suggests that N is 
not highly enriched. The various, general properties of 
RP916 are summarized in Tabled 

The original point of our work was to search for vari- 
ability in planetary nebulae. This search was enabled by 
using the automated pipeline developed for the Super- 
MACHO project, which makes use of a very powerful 
difference irriage t echnique (see Alard fc LuptonI 1998t 
Becker et al.l 120041 ). In brief, all of the images are re- 
sampled to the same geometry and sky subtracted, and 
the images are placed on the same photometric scale. 
The PSF of the template image is then matched to each 



target image via convolution, and subtracted. The re- 
sult of this subtraction on four selected nights is shown 
in Figure [HI where the difference has the value 0.0 every- 
where (apart from shot noise) except for features that 
are either fainter or brighter than those in the template 
image. Photometry is then performed at the positions 
of the lobes on the difference image for each observation. 
This technique has been shown to yield accurate pho- 
tometry even in crowded fields such as the LMC, and is 
extremely robust against non-photometric conditions. 

Light curves were determined for the eastern and 
western lobes of RP916, which are shown in Figured! 
Note that both lobes increased in brightness compared 
to the template image, but some time after the second 
year of the observing campaign the eastern lobe began 
to fade while the western lobe continued to brighten. 
By MJD - 53, 200 = 2004 July the western lobe began 
to fade, but more slowly than the eastern lobe. Close 
inspection of Figure [3] suggests that the orientation of 
the brightness enhancement in the western lobe changed 
and became elongated by MJD 53, 735 = 2005 Dec. 
The western lobe brightened by nearly 40% during the 
course of the observing campaign. 

4. Discussion 

RP916 is unusual in many respects, including its large 
physical size, strong IR emission, and optical variability. 
But the morphology and moderate ionization spectrum 
leave little doubt that it is a genuine planetary nebula, 
and the measured radial velocity establishes its location 
in the LMC. If the central star of RP916 has a binary 
companion that is very red and luminous it could con- 
tribute to the near-IR continuum luminosity, but it can- 
not account for the extreme colors in the mid-IR. Note, 
however, that since the epochs of the 2MASS and SAGE 
observations differ by 5.4 years it is possible that the ra- 
tio of the near- to far-IR luminosity could have changed 
somewhat during that time. We conclude that RP916 
contains a large amount of dust (which might be ex- 
pected given the pinched waist and the dark band that 
bisects the nebula). If RP916 is similar to other LMC 
PNe with significant dust, both molecular emission fea- 
tures and thermal continuum radiation contribute sub- 
stantially to the total IR flux, part icularly in the mid-IR 
(see, e.g.. lStanghenini et al.ll2007 ). 

The nebular variability is the most unusual feature 
of this object, and it is difficult to think of a cause that 
does not involve a binary star as the central engine. In- 
deed, the lack of an extreme N e nrichment, which is 
typical of bipolar PNe in the LMC (|Shaw et al.ll2006t ) is 
consistent with common-envelope evolution, where the 
expected conversion of C to N via hot-bottom burning 
is suppressed. It may be that the variation resulted 
from an outburst from a red companion. In this case 
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a change in optical luminosity might be manifested in 
the nebula as a light echo seen in scattered light from 
dust within the nebula. (Note that the phase lag of 
~ 300 — 500 d between the peak emission of the two 
lobes in Figure 2] would imply a significant inclination 
out of the plane of the sky.) However, no Mira-like peri- 
odicity is detected over the 1450 d observing campaign. 
It seems unlikely that the variation is from a global 
change in ionizing photons, as the inferred H recom- 
bination timescale is ~ 300 yr if the derived density is 
representative of the whole nebula. In any case, a light 
echo cannot explain the radial structures that are evi- 
dent in the detailed HST images. The variability might 
on the other hand result from a precessing jet of ma- 
terial from the central source that interacts with (i.e., 
shocks) the surrounding nebula. Other PNe that show a 
similar photometric variability, resulting from a Symbi- 
otic com panion, are the Galactic bipolar P N He2-104 



tCorra di et al. 

mm 



20011) and LMC-SMP 83 (|Peiia et all 



Dovle et al. (|2000l) interpreted the variability 



within M2-9 as an interaction of a jet from the cen- 
tral source interacting with the surrounding PN, which 
has over time continued to shape the nebular morphol- 
ogy. While this mechanism is an interesting possibility 
for RP916, confirmation would require long-slit spec- 
troscopy to measure the velocity field, high-resolution 
(~ O'.'l or better) imaging in Ha to reveal the time- 
dependent morphology of the variability, imaging po- 
larimetry to measure the extent of any scattered light 
by dust, and infrared spectra to understand the nature 
of the IR luminosity and (possibly) detect a red, lumi- 
nous companion star. 

The frequency of nebular variability generally among 
PNe is not known, but RP916 may be the most ex- 
treme example of how PNe morphology can be re-shaped 
even at very advanced stages of evolution. Photometric 
monitoring at high spatial resolution of a sizable sample 
of bipolar PNe in the LMC is crucial for determining 
what fraction of PNe experience the unusual variability 
of RP916. Perhaps the observing campaigns of the next 
generation of survey telescopes such as that for LSST, 
which will have the necessary depth, sky coverage, and 
cadence, can help to resolve this broader question. 
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Table 1 
IR Photometry 





SAGE'' 






2MASS'' 


Band 


Magnitude 


Flux (mjy) 


Band 


Magnitude 


S.dfim 


12.46 ±0.15 


2.9 ±0.4 


J 


16.32 ±0.18 


4.5/im 






H 


15.07 ±0.12 


5.8^m 


10.35 ±0.04 


8.5 ±0.3 


K 


13.88 ±0.08 


S.OAtm 


8.58 ± 0.03 


23.2 ± 0.7 






2Afj.m 


3.34 ±0.01 


340. ± 3. 







^SAGE project catalog entry J053059.45-683541.2 



llMeixner et al.ll2006t) 

^2 MASS catalog entry 05305947-6835413 l|Skrutakie et al.l 
l2006tl 



Table 2 
Basic Properties of RP916 



Property 


Value 


Reference 


R.A. (J2000) 


5*' 30™ 59? 481 


1 


Dec (J2000) 


-68° 35' 41'.'29 


1 


Morphology 


"Butterfly" bipolar 


1, 2 


Angular Size 


7'.'0 X 3'.'0 


1 


Physical Size 


1.72 X 0.74 pc 


1 


Density 


400 cm" 3 


3 


yHelio 


280 km s-l 


2 


c 


0.9 


3 



Ref eren ces. — 1. this work; 
1 I2OO6D : 3. lReidl tOOt i 



2. iReid &: Parked 
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Fig. l.~ Images of RP916 in the SuperMACHO VR 
bandpass (left) with square-root intensity scale, and in the 
AAO/UKST I3a + R continuum (right) with log intensity 
scale. Images are 30" x 30" with north up and east to the 
left. The image resolution in VR is ~ 0"8, and in Hq + R 
is ~4". The nebula is marked in the Ha + R image with 
a green, horizontal bar; pure emission lines are rendered as 
blue, and continuum sources are rendered in pink. 




Fig. 3. — Difference between the template image of RP916 
taken on 2002 Jan. 16 = MJD 52290 (Fig.© and those taken 
on 2002 Dec. 14 = MJD 52622 (upper left), 2003 Dec. 19 = 
MJD 52992 (upper right), 2004 Dec. 13 = MJD 53352 (lower 
left) and on 2005 Dec. 31 = MJD 53735 (lower nght). Image 
intensity scale is linear, with zero corresponding to medium 
grey. Image size and orientation as in Fig.[T] Saturated areas 
in the input images are set to zero. 




Fig. 2.— HST image of RP916 in the F814W bandpass, 
with log intensity scale. Image is 20" x 20" with north up 
and east to the left. 



Fig. 4. — Light curves for eastern (triangles) and western 
(circles) lobes of RP916. Abscissa shows the SuperMACHO 
VR bandpass. The epoch of the template image is indicated 
(solid line), as is the epoch of the optical spectrum (dashed 
line) and the SAGE mid-IR campaign (dotted lines). 
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